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PARAMETRIC INVESTIGATION OF LIQUID HYDROGEN 
TANK PRESSURIZATION DURING OUTFLOW 
by David A. Mandell and William H. Roudebush 
Lewis Research Center 

SUMMARY 

The factors influencing the amount of gas required to pressurize a cylindrical tank 
of liquid hydrogen during outflow are investigated analytically based on a simplified 
model of the tank-flow and heat-transfer problem. The investigation begins with an es- 
tablished set of differential equations describing the problem. These equations and as- 
sociated initial and boundary conditions are put in terms of dimensionless variables. 

The ratio of the mass of gas required for pressurization to a standard mass is shown to 
depend on only a few parameters. By computer solution of the equations over a wide 
range of values of the parameters it is found that only two are of principal importance. 
These two dimensionless parameters have the form of modified Stanton numbers, one 
associated with the gas and one with the tank wall. A figure is presented from which the 
pressurant mass ratio can be found for a wide range of values of the modified gas and 
wall Stanton numbers. A comparison with experimental results is made. 


INTRODUCTION 

Liquid hydrogen is now established as a rocket fuel. However, the design of new 
hydrogen rocket systems will require a better understanding of the thermal processes 
occurring inside the hydrogen tank. The particular problem treated here is the predic- 
tion of the amount of pressurizing gas that must be added to the tank to maintain a de- 
sired pressure during outflow of the fuel. 

A one-dimensional analysis of the nonsteady heat-transfer and flow problem for a 
cylindrical tank during outflow was given in reference 1. A numerical method was pre- 
sented there for calculating the temperature distribution in the tank, from which the 
amount of pressurizing gas that must be supplied during outflow can be obtained. Cal- 
culated results in reference 1 were compared with experimental data. The agreement 



was sufficiently good to suggest the use of the method for an analytical investigation of 
the various parameters influencing the pressurant requirement. 

Such an investigation is described in this report. It is shown that the important 
factors entering the pressurization problem can be combined into a small number of 
parameters. These parameters are then varied over a wide range and their effect on 
mass requirements is obtained. 

The analytical results presented herein extend far beyond the range where the 
theory of reference 1 has been checked against experimental data. Furthermore, some 
complicating factors such as liquid sloshing and tank dome heating are ignored. Never- 
theless, it is hoped that the results will give at least a reliable estimate of the effects of 
the most important pressurization parameters for cylindrical tanks and help in prelimi- 
nary design. 


ANALYSIS 

Assumptions 

The following parametric analysis of the pressurization problem is based on the 
equations of reference 1. In that reference a simple flow model was constructed based 
on the following assumptions: 

(1) The ullage gas is nonviscous. 

(2) The ullage gas velocity is everywhere parallel to the tank axis and does not vary 
radially or circumferentially. 

(3) The tank pressure does not vary spatially. 

(4) The ullage gas temperature does not vary radially or circumferentially. 

(5) The tank wall temperature does not vary radially or circumferentially. 

(6) No heat is transferred axially in either the gas or the wall. 

(7) No condensation or evaporation occurs. 

(8) No heat is transferred at the liquid interface or at the top of the tank. 

With these assumptions the problem reduces to a one-dimensional nonsteady flow prob- 
lem with heat addition. The governing fluid dynamic and heat-transfer equations are 
given in reference 1. 

Some further assumptions are made in the present report. These added restric- 
tions allow the differential equations and the initial and boundary conditions to be put 
into a simple form that clearly exhibits the significant dimensionless groups. Later in 
the report an examination is made of the effect of these assumptions: 

(9) The ullage gas is a perfect gas with constant specific heat. 
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(10) The inlet gas temperature, the tank pressure, and the outflow rate are constant. 

(11) The gas-to-wall heat -transfer coefficient is constant. 

(12) The gas and wall temperatures at the start of outflow are equal and vary 
linearly in the direction of the tank axis from saturation temperature at the liquid inter- 
face to one-half the sum of the saturation temperature and the inlet gas temperature at 
the top of the tank. 

( 13) The gas and wall temperatures at the liquid interface are constant and equal 
throughout the outflow period. 


Differential Equations and Initial and Boundary Conditions 

A schematic drawing of a cylindrical tank is shown in figure 1. With the added as- 
sumptions (9) to (13) the equations of reference 1 for such a tank reduce to 
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Figure 1. - Schematic of cylindrical liquid- 
propellant tank during outflow. 


All symbols are defined in the appendix. The detailed 
development of these equations and a discussion of the 
underlying assumptions (1) to (8) are given in refer- 
ence 1, along with a comparison of calculated and ex- 
perimental results. 

Let L(t) denote the x-location of the liquid sur- 
face at time t (see fig. 1). The velocity u L of the 
surface is constant since the outflow rate has been as- 
sumed constant. Therefore, 

L(t) = L 0 + u L t (4) 

where L Q s L(0) is the axial length of the ullage at the 
start of outflow. 

From assumption (12) the initial conditions can be 
written as 
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T = T W = TL + Tg + — ( TL ~ Tg ) t = 0 0^x£L o (5) 

W 2 L 0 \ 2 / 

This choice of initial temperatures is arbitrary. In an actual example the initial gas and 
wall temperatures need not be equal and need not relate in any way to T . The effect of 
the choice will be examined later in the report. 

The gas temperature at the top of the tank during outflow is 

T = T = constant x = 0 0 < t ^ t f (6) 

o 

The gas and wall temperatures at the liquid surface are 

T = T w = Tj^ = constant x = L Q + u^t 0 < t S t^ (7) 

and the velocity of the gas is 

u = u-^ = constant x = L Q + u^t 0 < t 2 t^ (8) 

Equations (1) to (8) completely determine the problem. 


Dimensionless Equations 


The preceding equations will now be made dimensionless by introducing the following 
dimensionless variables: 
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Since the velocity of the liquid surface is constant, 


L f" L o _u L t f (10) 

Using equations (9) and (10) in equations (1) to (8) yields the following equations: 
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When the perfect gas equation of state and equation (10) are used, the group of quan- 
tities appearing in equation (11) can be written as a modified Stanton number: 


hRt,T 
f g 

rMPCp 


ht. 


rp c„ rp c u T 
H S P g P L 


= St 


g 


(19) 


The modified gas Stanton number St^ is a true Stanton number modified by division by 
the dimensionless radius ? . 

In connection with equation (12) a second modified Stanton number related to the wall 
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and called the modified wall Stanton number is defined as 


htr j_ 

St_ = £ = * (20) 

^w^ 3 w C w w C w u L 

Although St w has the form of a Stanton number, properties of both wall and fluid are in- 
volved. The use of such a dimensionless group is unusual. An alternative is to replace 
St w by the parameter 


St g Vw c w 


which is one-half the ratio of the heat capacity of the wall to the heat capacity of the gas. 
However, St w has been retained in this analysis since it appeared naturally in the dif- 
ferential equation. 

An examination of equations (11) to (20) reveals that the set of dimensionless dif- 
ferential equations and the initial and boundary conditions depend only on the four dimen- 
sionless quantities St , St , L , and TV • With the assumptions previously made for 

a . o ™ u •* -J 

this analysis, St^, L q , and are constant for a particular problem. If, furthermore, 
c w can be considered constant, the following result is obtained: Every solution of the 
pressurization problem in terms of the dimensionless variables (eq. (9)) is completely 

a 

determined by the specification of the four dimensionless constants St^, St w , L q , 
and T^. 

In dealing with liquid hydrogen problems the wall specific heat variation may be 
significant. In that case the parameter St w will vary, also. In order to maintain St w 
as a characteristic constant of the problem, it is convenient to introduce some reference 
value of specific heat U w constant for a given problem. To do this equation (20) is 
modified in the following way: 
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Figure 2. - Integrated average wall specific heats used 
to calculate waif Stanton numbers for stainless- 
steel and aluminum tanks. 


In equation (21) the parameter St is constant for a given problem, and the variation of 

W A 

wall specific heat is accounted for in the new dimensionless variable c . The reference 

w 

specific heat c can be chosen arbitrarily; for example, a value at room temperature 
w 

could be used. However, the use of an appropriate average value would be expected to 
reduce the effect of the variable c . In this report the following value is used: 



A plot of c against T is given in figure 2 for stainless steel and aluminum. 

W o 

When equation (21) is used, the problem is seen to depend on the four dimensionless 

^ A A 

constants Stg, St w , L q , and and the dimensionless specific heat variation c w . 

The specific heat variation c w can be expected to vary differently with different 
tank wall materials and thereby produce different results. Therefore, attention will be 
confined largely to a single tank wall material, stainless steel. Even when this is done, 
however, a difficulty remains. The variation c w is a function of the dimensional var- 
iable T as well as the dimensionless variable T. Thus, the following result can be 
formulated: For a given wall material every solution of the pressurization problem in 
terms of the dimensionless variables of equation (9) is completely determined by the 
specification of four dimensionless constants Stg, St w , L q , and T^ and the one- 

dimensional constant T . 

g 
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Mass Ratio 


A A. A 

A solution of the dinaensionless pressurization equations gives T, T w , and u as 
functions of x and t. When the assumption of a perfect gas is used, the total mass of 
pressurizing gas added during outflow is given by 


“ = I" 


* 7rr^p(x, tj)dx - / ° 


7rr p(x, 0)dx 


nr^MP 

R 




7rr^MP(L^ - L Q ) 


RT 


g 



(23) 


If the ideal pressurant mass nij is defined to be the mass required to pressurize the 
tank if no heat transfer occurs, then 


m^MP(L f - L ) 

m i = ^^f - L o )p g = ^rT- 

g 

Combining equations (23) and (24) gives 



(24) 


(25) 


The dimensionless temperature variation T(x, t) appearing in equation (25), and thereby 
the mass ratio m/nr , is determined by the solution of the dimensionless equations. 
This mass ratio is sometimes referred to as a collapse factor in pressurization litera- 
ture. As just shown it arises naturally from the choice of dimensionless variables in 
equation (9). 

The mass ratio but not the mass itself is determined by the solution of the dimen- 
sionless equations. Determining the mass from the mass ratio requires also the dis- 
placed liquid volume 7 rr^(L f - L Q ), the tank pressure, and the inlet gas temperature 
(see eq. (24)). 
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Mass ratio, m/m^ 




Figure 3. - Reference Stanton number map showing values of pressurant mass ratioor range of gas and w 

modified Stanton numbers for cylindrical stainless-steel tanks. Inlet gas temperature. 500 R. liquid sur 

surface temperature, 37° R; dimensionless initial ullage height, 0.0526. 
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Wall Stanton number, St w 
(a) Gas Stanton number, 5. 



Gas Stanton number, Stg 
(b)Wall Stanton number, 2.5. 


Figure 4. - Effect on mass ratio of changing the inlet gas temperature. Dimensionless interface temperature, 

0.074; dimensionless initial ullage height, 0.0526. 

RESULTS 

As indicated in the previous section, the mass ratio m/im is completely deter- 
mined by choosing a tank wall material and fixing values of the five constants St^, St^, 
L Q , T^, and T . In the present section some quantitative results are presented for 
stainless -steel tanks. 

First, with the computer solution described in reference 1, values of gas tempera- 
ture distribution and thereby values of m/m. are computed for a wide range of values of 
Stg and St w by using fixed values of L q , T-^, and T^. Then with these results as a 
standard for comparison, additional calculations are made to examine the effect of T , 
T^, and L q . Finally, some of the assumptions made to obtain the simple dimensionless 
equations are relaxed, and calculations are made to examine their effect. 
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(a) Gas Stanton number, 5. 



Figure 5. - Effect on mass ratio of changing the dimensionless interface temperature. Inlet gas temperature, 
500° R; dimensionless initial ullage height, 0.0526. 


Reference Stanton Number Map 

The inlet gas temperature T is taken to be 500° R. The temperature Tj is 

chosen to be 0. 074, corresponding to a liquid surface temperature of 37 R. Selection 

of a 5-percent initial ullage volume gives a value of L q = 0. 0526. The effect of this 

choice of values for T , T t , and L will be examined later. With these values fixed, 

g* Lr O 7 

the mass ratio becomes a function only of the modified Stanton numbers. For each pair 
of values of St^ and St w a computer solution can be carried out to obtain m/ rm. 

The computed mass ratios for a wide range of Stanton numbers are shown in figure 3. 
This figure will be referred to as the reference Stanton number map. The map is in- 
tended to encompass existing experimental tanks and flight tanks that now exist or are 
likely to be developed. In the next section it is shown that the effects of T , T L , and 
are relatively small so that figure 3 provides a rapid estimate of mass ratio for a 
wide variety of cylindrical tanks. 
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(a) Gas Stanton number, 5 



Gas Stanton number, St g 
(b)Wall Stanton number, 2.5. 

Figure 6. - Effect on mass ratio of changing the dimensionless initial ullage height. Inlet gas temperature, 
500° R; dimensionless interface temperature, 0.074. 


Effect of Parameters Tg, T L , and L 0 


When L and T t are held fixed, the mass ratio depends only on T for fixed 
values of the gas and wall Stanton numbers. As previously explained, this dependence 
on T enters through the dependence of the wall specific heat ratio c w on T . It 
might be expected therefore that the dependence is small. An indication of the effect of 


changing T from the value of 500° R used in the reference Stanton number map to 

300° and 700° R is given in figure 4. At a value of St = 5 the maximum deviation from 

& 

the reference map values of mass ratio is 7. 1 percent. At higher values of St g the 


deviation increases as shown. 


The effect of changing while holding L q and fixed at reference map values 
is shown in figure 5. To change T L the liquid surface temperature T L is changed 
from 37° to 60° R. The effect appears to be negligible except at large values of St g . 

A similar result is true for L Q . Keeping T L and T fixed at reference map 
values and changing L Q gives the result shown in figure 6. The difference in mass 
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(a) Gas Stanton number, 5. 



Gas Stanton number, Stg 
(b) Wall Stanton number, 2.5. 


Figure 7. - Effect on mass ratio of allowing gas specific heat to vary with temperature. Inlet gas temperature, 
500° R; dimensionless interface temperature, 0.074; dimensionless initial ullage height, 0.0526. 


ratio brought about by changing the initial ullage volume from 5 to 20 percent is less 
than 7 percent over the range of Stanton numbers shown in the figure. When the initial 
ullage volume is increased to 50 percent, the maximum difference in mass ratio in- 
creases to about 14 percent. 

These calculated results indicate that St^ and St w are the primary factors affect- 
ing the mass ratio, at least for moderately small values of initial ullage. 


Effect of Certain Assumptions 

The effect of assuming gas specific heat to be constant was tested by computing mass 
ratios with a variable specific heat. In order to make a comparison with the values ob- 
tained for constant specific heat, it is necessary to assign a value of Stg to the cases 
where specific heat is allowed to vary. This is done by computing a mean c defined 

by 

1 P« 

c D = ± I c (T)dT 

p v 37 L 
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Wall Stanton number, St w 
(a) Gas Stanton number, 5 . 



Gas Stanton number, St g 
(b) Wail Stanton number, 2.5. 


Figure 8. - Effect on mass ratio of holding wall specific heat constant Inlet gas temperature, 500° R; dimension- 
less interface temperature, 0.074; dimensionless initial ullage height, 0.0526. 

and by using it to determine St for each case. With this done the constant and varying 

o 

specific heat cases can be compared at equal Stanton numbers. The results, shown in 
figure 7, indicate a negligible difference in mass ratios. 

A variable wall specific heat was used to obtain the reference map (fig. 3). Fig- 
ure 8 shows the effect of using instead a constant value of specific heat taken to be the 

integrated average value of figure 2. At St g = 5 a maximum difference in mass ratio 
of 18. 5 percent occurs at St w = 1. Even larger differences occur at larger values of 
St . It was because of these large differences that the wall specific heat was allowed to 

IS 

vary for the reference map. 

It is interesting, however, that changing the wall material from stainless steel to 
aluminum does not affect the mass ratio, as shown in figure 9. The reason for this is 
found in the near similarity of the variations c w with T w for the two metals (fig. 10). 

Another assumption needed for the simplification of the dimensionless problem was 
that the wall and gas temperatures are equal and vary linearly with distance at the start 
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(a) Gas Stanton number, 5. 


1 1 1 

.01 .02 

(b) Wall Stanton number, 2.5. 

Figure 9. - Effect on mass ratio of changing tank wail material from stainless steel to aluminum. Dimensionless 
interface temperature, 0.074; dimensionless initial ullage height, 0.0526. 




0 .2 .4 .6 .8 ^ 10 

Dimensionless wall temperature, T w 


Figure 10. - Dimensionless wall specific heat 
as function of dimensionless wall tempera- 
ture for stainless steel and aluminum. 

Inlet gas temperature, 500° R. 
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Map value, T = T w 2 ^ at t = °- x - *-o 


31 .02 
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-o-- T W = T L , T=(T L -Tg) ^ +T g at t = 0, 0 < x < 1^ 
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Wall Stanton number, St w 


(a-1) Gas Stanton number, 5. 


.02 



(a-2) Wall Stanton number, 2.5. 


(a) Dimensionless initial ullage height, 0.0526. 



Wall Stanton number, St w 
(b-1) Gas Stanton number, 5. 



Gas Stanton number, St g 
(b-2) Wall Stanton number, 2.5. 

(b) Dimensionless initial ullage height, 0.25. 

Figure 11. - Effect on mass ratio of changing the assumed initial wall and gas temperatures. Dimensionless 
interfacetemperature, 0.074. 




(b)Wall Stanton number, 2.5. 

Figure 12. - Effect on mass ratio of considering a real gas. inlet gas temperature, 500° R; dimensionless inter- 
face temperature, 0.074; dimensionless initial ullage height, 0.0526. 

of outflow. In particular, they were assumed to vary from at the liquid surface to 
<T l + T )/ 2 at the top of the tank. Since the initial temperatures may not be easily 
computable in a particular case, this assumption assumes considerable importance. 

To determine how critical the initial temperature distribution is, the following 
change in initial values was made. The initial wall temperature was taken to be 37° R at 
all points, and the gas temperature was varied linearly from 37° R at the liquid surface 
to 500° R at the top of the tank. The effect that this severe change in initial conditions 
had on the mass ratio is shown in figure 11. For initial ullage volumes of 5 percent 
(corresponding to values of ~L Q j - L Q ) = 0. 0526) the deviation from reference map 
values of mass ratio is negligible. At initial ullages of 20 percent the deviation is more 
significant, increasing to about 12 percent at St^ = 5 and St w = 0. 1. 

The effect of assuming a perfect gas was also tested by using a real gas compress- 
ibility factor Z(T, P) for hydrogen. The details of the equations needed in this case are 
given in reference 1. The results shown in figure 12 indicate that the difference in mass 
ratio brought about by using a real gas in the calculations is negligible. 

It is difficult to evaluate the validity of the assumption of constant heat-transfer co- 
efficient since the correct variation is not known for comparison. Reference 1 shows, 
however, that good agreement (within about 5 percent on the average) with experimental 
mass ratios can be obtained by using a constant value of h in the calculations. The 
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(b)Wall Stanton number, 2.5. 

Figure 13. - Effect on mass ratio of transient tank pressure variations. Inlet gas temperature 500° R; dimension- 
less interface temperature, 0.074; dimensionless initial ullage height, 0.0526. 


constant value used in each case in reference 1 was an average of experimentally mea- 
sured local values of h. 

The assumption of constant values of tank pressure, inlet gas temperature, and 
outflow rate will be violated to some extent in all real situations. In many cases, 
however, the deviation from constant values of these quantities will occur primarily 
during a starting transient. Figures 13, 14, and 15 show the effect on mass ratios of 
starting transients in tank pressure, inlet gas temperature, and outflow rate, respec- 
tively. Transients of 10 and 50 percent of total outflow time were considered. The 
magnitude of the transients are shown in each case in the respective figure. 

The pressure transient has the most pronounced effect. The deviation, for all 
values of St^ and St w , is about 25 percent for the 50-percent transient pressure varia- 
tion and about 10 percent for the 10-percent transient. 

For the 50-percent transient gas temperature variation, the deviation from the map 
value is approximately 14 percent at St^ = 5 and St^ = 0. 01. This deviation decreases 
to 0 percent at St^ = 5 and St w = 100. For St w = 2. 5, the deviation is approximately 
14 percent at St_ = 100 and decreases to 1 percent at St = 1. The deviation is less 

5 & 
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(a) Gas Stanton number, 5. 



Gas Stanton number, St g 
(b)Wall Stanton number, 2.5. 


Figure 14. - Effect of transient inlet gas temperature variations. Inlet gas temperature, 500° R; dimensionless 
interface temperature, 0.074; dimensionless initial ullage height, 0.0526. 



(a) Gas Stanton number, 5. 



Gas Stanton number, St g 
(b) Wall Stanton number, 2.5. 

Figure 15. - Effect on mass ratio of transient liquid outflow variations. Inlet aas temperature, 500° R; dimension- 
less gas ratio, 0.074. 
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(b) Wall Stanton number, 2. 5. 

Figure 16. - Final dimensionless gas temperature variation at end of outflow. Dimension- 
less initial ullage height, 0.0526; dimensionless interface temperature, 0.074; dimen- 
sionless time from start of outflow, 1. 


than 4 percent for the 10-percent transient. 

The 50-percent transient liquid outflow variation deviates from the map a maximum 
of 8 percent. The 10-percent transient deviates about 2 percent. 


Dimensionless Temperature Distributions 

A 

In addition to the propellant mass ratios, the dimensionless temperatures Tg and 
are given as functions of x and t in the computer solutions. Figures 16 and 17 
show the variety of shapes that can occur in the final gas and wall temperature profiles, 
respectively, obtained at various points of the Stanton number map. The assumption 
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(b)Wall Stanton number, 2.5. 


Figure 17. - Final dimensionless wall temperature at end of outflow. Dimensionless ini- 
tial ullage height, 0.0526; dimensionless interface temperature, 0.074; dimensionless 
time from start of outflow, 1. 


sometimes used in simple analyses of the pressurization problem, that a linear variation 
of temperature can be used as a good approximation to the actual temperature distribu- 
tion, is evidently not valid over the entire range of Stanton numbers. 


COMPARISON OF PREDICTED AND EXPERIMENTAL 
VALUES OF MASS RATIO 


In obtaining the reference Stanton number map (fig. 3, p. 7) many simplifications 
were made. The foregoing results indicate, however, that the most important effects on 
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TABLE I. - VALUES USED FOR COMPUTATION OF STANTON 


NUMBERS FOR EXPERIMENTS OF REFERENCE 2 
[Wall thickness, 0.026 ft; wall density, 500 1b/cuft; tank radius, 1.13 ft.] 


Run 

Average wall 
specific heat, 

c w’ 

Btu 

Gas density, 

P g’ 

lb 

cu ft 

Gas specific 
heat, 

V 

Btu 

Outflow 

time, 

h* 

sec 

Heat -transfer 
coefficient, 
hx3600, 
Btu 

(lb)(°R) 

(lb)(°R) 

(hr)(ft2)(°R) 

1 

0.0730 

0.0573 

2. 

98 

350 

13. 75 

2 

. 0725 

.0584 



93 

12.25 

3 

. 0730 

.0204 



284 

7. 09 

4 

.0735 

. 0206 



101 

6.67 

5 

.0423 

. 1110 



95 

11.34 

6 

.0800 

.0120 



88 

5. 13 

7 

.0738 

. 1120 

1. 24 

355 

12.31 

8 

.0738 

. 1120 



90 

11. 15 

9 

.0316 

.2790 



100 

10.45 

10 

.0792 

. 0240 



309 

5. 25 


the pressurant mass ratio are probably included in the map. To test this experimen- 
tally, the liquid hydrogen outflow experiments that were used for comparison in refer- 
ence 1 will be considered here. These experiments (reported in refs. 2 and 3) cover a 
range of tank pressure, outflow rate, and inlet gas temperature, and include, in some 
instances, significant transients in the inlet gas temperature. 


Reference 2 Data 

In reference 2 the authors report some of the results of a systematic series of 
liquid hydrogen expulsion experiments with a cylindrical tank. The tank used was 
27 inches in diameter and 89 inches in overall length with dished head ends. A gas dif- 
fuser was used at the inlet. The tank was made of 5/16-inch 304 stainless-steel plate 
and was vacuum jacketed. The internal instrumentation, described in detail in refer- 
ence 2, provided a relatively significant heat sink in some of the tests. 

Ten experimental runs (not all of which were published in ref. 2) were used in refer- 
ence 1 to compare calculated and experimental mass requirements. These rims covered 
a wide range of values of temperature, pressure, and outflow rate. Four of the runs 
used helium as the pressurizing gas. The detailed data are given in reference 1. Table I 
shows the data needed for the calculation of modified gas and wall Stanton numbers for 
each case. The values of heat-transfer coefficient are average values determined from 
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TABLE n. - COMPARISON OF CALCULATED AND 
EXPERIMENTAL MASS RATIOS FOR 
EXPERIMENTS OF REFERENCE 2 


each experiment. The wall specific 
heat is the integrated average value 
indicated in figure 2 (p. 7). 

When the values in table I were 


Run 

Dimension- 
less gas 
Stanton 
number, 
St 

g 

Dimension- 
less wall 
Stanton 
number. 

Mass rat 
Figure 3 

lo, m/m. 

Experi- 

mental 

Percent 

difference 

1 

6.93 

1.41 

2. 58 

2.97 

-13. 1 

2 

1. 61 

.34 

1. 77 

2. 12 

-16.5 

3 

8. 14 

.59 

3.09 

3.88 

-20.3 

4 

2. 70 

. 20 

2. 20 

2. 66 

-17.3 

5 

00 

00 

. 54 

1.47 

1. 53 

-3.9 

6 

3. 10 

. 12 

2. 37 

3.08 

-23.0 

7 

7. 73 

1. 27 

2. 73 

3. 24 

-15.7 

8 

1. 78 

. 29 

1. 86 

2. 14 

-13. 1 

9 

. 74 

.71 

1. 38 

1.41 

-2. 1 

10 

13. 39 

.44 

3.92 

5. 10 

-23. 2 


used, the modified gas and wall 
Stanton numbers were computed 
from equations (19) and (21), and the 
values of mass ratio were then 
found from figure 3. Table II lists 
the computed values of modified 
Stanton numbers and the resulting 
values of mass ratio. Also listed 
are the experimental values of mass 
ratio in each case and the percent 
difference between the calculated 
and experimented values. 

The agreement shown in table II 
is not as good as might be expected 
from the parametric analysis. The 


explanation for this probably lies in 
the large heat sink provided by the internal instrumentation in the experiments of refer- 
ence 2. The effect of the heat sink on mass ratio, which is not accounted for in figure 3, 


was shown to be significant in the calculations of reference 1. In that reference the ef- 
fect of the heat sink on the mass requirement was determined from an estimated value of 
heat flow to the internal hardware. If the experimental values of m/im given in table II 
are adjusted by these calculated amounts, the results are those given in table HI. With 
the adjustment for the heat flow to the internal hardware, the agreement of the calculated 
and experimental values is much better. The average absolute value of the difference 
for the 10 cases is 6. 2 percent. It is clear here, as it was in reference 1, that neglect- 
ing internal heat flow can cause significant error in the calculated results. 


Reference 3 Data 

Reference 3 reports the results of liquid hydrogen expulsion experiments using a 
40-inch-diameter cylindrical test tank 100 inches in overall length. The test tank was 
0. 090-inch-thick stainless steel and was enclosed in a 60-inch-diameter vacuum-tight 
carbon- steel tank. A gas diffuser was in the top and an antivortex baffle was in the 
bottom. Perforated conical slosh baffles were located at various axial distances. The 
heat sink effect of the internal hardware could not be estimated from the data reported 
in reference 3. 
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TABLE HI. - COMPARISON OF 
CALCULATED AND ADJUSTED 
EXPERIMENTAL MASS RATIOS 


TABLE IV. - VALUES USED FOR COMPUTATION 
OF STANTON NUMBERS FOR EXPERIMENTS 
OF REFERENCE 3 


FOR EXPERIMENTS OF 
REFERENCE 2 


Run 

Mass ratio, m/m^ 

Percent 


Figure 3 

Adjusted 

experi- 

mental 

difference 

1 

2. 58 

2. 76 

-6. 5 

2 

1. 77 

1.84 

-3. 8 

3 

3.09 

3. 31 

-6. 7 

4 

2. 20 

2. 25 

-2. 2 

5 

1.47 

1. 35 

8.9 

6 

2. 37 

2. 56 

-7. 4 

7 

2. 73 

2, 86 

-4. 5 

8 

1.86 

1.93 

-3. 6 

9 

1. 38 

1. 25 

10.4 

10 

3.92 

4. 25 

-7. 8 


[Wall thickness, 0.0075 ft; wall density, 500 lb/cu ft; 
tank radius, 1. 67 ft.] 


Run 

Average 

Gas 

Gas 

Out- 

Heat- 


wall 

density, 

specific 

flow 

transfer 


specific 


heat, 

time, 

coefficient, 


heat, 

lb 

V 

tf, 

hx3600, 


c w’ 

Btu 

cu ft 

Btu 

sec 

Btu 


(lb)(°R) 


(hr)(ft 2 )(°R) 


(lb)(°R) 




1 

0. 0475 

0.0283 

2. 67 

89 

11. 5 

2 

.0730 

. 0171 

2.98 

103 

a 12. 0 

3 

.0475 

.0289 

2. 67 

120 

11.3 

4 



. 0289 

2. 67 

87 

12.0 

5 



. 0566 

1. 24 

99 

12. 1 

6 



.0280 

2. 67 

111 

11. 8 

7 



. 0286 

2. 67 

97 

11. 7 

8 



. 0286 

2. 67 

105 

13. 9 


a Estimated value - not given in reference 3. 


TABLE V. - COMPARISON OF 
CALCULATED AND EXPERI- 
MENTAL MASS RATIOS 
FOR EXPERIMENTS OF 
REFERENCE 3 


Run 

Mass ratio, m/m. 

Percent 


Figure 3 

Experi- 

difference 



mental 


1 

1.72 

1. 54 

11. 7 

2 

2. 14 

2. 16 

-.9 

3 

1. 79 

1. 61 

11. 2 

4 

1. 71 

1. 68 

1.8 

5 

1. 81 

1. 69 

7. 1 

6 

1. 80 

1. 71 

5. 3 

7 

1.75 

1. 59 

10. 0 

8 

1.83 

1. 69 

8. 3 


Nine tests were reported in reference 3 for which the 
system vacuum was maintained. Of these nine tests (all 
of which were used in the comparisons of calculated and 
experimental results in ref. 1) only eight are considered 
here. It was necessary to omit the one that used hydro- 
gen to pressurize the tank but had helium in the ullage 
initially. Such an arrangement is not covered in the re- 
sults of figure 3. The cases used cover two values of in- 
let gas temperature and a range of initial ullage volumes 
up to about 20 percent. The outflow time and tank pres- 
sure vary slightly for the eight cases. Helium was used 
as a pressurizing gas in one case. Sloshing was induced 
in all but one case. Complete data for the calculation are 
given in reference 1. 

Table IV shows the data used for the calculation of 
the modified Stanton numbers for the eight experiments. 
The values of h are experimental average values as they 
were in table I. 

When the values in table IV were used, the modified 
Stanton numbers were computed as before and the mass 
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TABLE VI. - COMPARISON OF MASS RATIOS 
FOR COMPUTED VALUES OF AVERAGE 
HEAT-TRANSFER COEFFICIENT 


Run 

He at -transfer 

Mass ratio, m/m i 

Percent 


coefficient, 

Figure 3 

Experi- 

difference 


equation (26) 


mental 



Reference 1 data 


1 

16.4 

2. 75 

2. 76 

-0.4 

2 

16. 5 

1.90 

1.84 

3.3 

3 

00 

co 

CJl 

3. 25 

3.31 

-1.8 

4 

8.40 

2. 35 

2.25 

4.4 

5 

18. 5 

1. 60 

1.35 

18. 5 

6 

6. 23 

2.45 

2. 56 

-4.3 

7 

12. 6 

2. 77 

2.86 

-3. 1 

8 

12. 6 

1.91 

1.93 

-1.0 

9 

17. 7 

1.45 

1.25 

16.0 

10 

4. 78 

3.80 

4. 25 

-10.4 


Reference 2 data 


1 

7.9 

1. 60 

1. 54 

3.9 

2 

7.4 

1.95 

2. 16 

-9.7 

3 

8. 1 

1.70 

1.61 

5.6 

4 

8. 1 

1. 55 

1. 68 

-7.7 

5 

6.9 

1. 65 

1. 69 

-2.4 

6 

7. 7 

1. 65 

1.71 

-3. 5 

7 

8.0 

1. 62 

1.59 

1.9 

8 

7.9 

1. 68 

1.69 

-. 6 


ratio found from figure 3 . Table V gives the 
resulting Stanton numbers and mass ratios 
and the comparable experimental values. The 
average absolute difference between calcu- 
lated and experimental values is 7. 0 percent. 

For data of references 2 and 3, there- 
fore, the agreement is quite good when proper 
allowance is made for internal heat sinks 
where these are important. This result bears 
out the implication of the parametric analysis 
that the reference Stanton number map con- 
tains the important factors influencing the 
mass ratio. It appears therefore that fig- 
ure 3 can be used to obtain rapid approximate 
predictions of mass ratio providing a suitable 
value of heat-transfer coefficient can be ob- 
tained. 

Heat -Transfer Coefficient 


The values of heat-transfer coefficient h 
used in determining Stanton numbers in the 
previous section were experimental values. 
For design use of figure 3 to obtain values of 
mass ratio, h must be estimated by some 
to use the free convection formula (ref. 4): 


other means. A simple way to do this is 


h = 0. 13 - (Gr X Pr)°- 333 

£ 


(26) 


where the properties are evaluated at the characteristic temperature Tg. The values of 
h obtained using equation (26) for the 18 experiments treated in the previous section are 
shown in table VI along with the experimental average values. Also shown in this table 
are the values of m/nu obtained from figure 3 by using the reference Stanton numbers 
calculated from the analytical values of h. The experimental values of m/nn are re- 
peated from tables IH and V and the final column gives the percent difference between 
calculated and experimental values of mass ratio. The average absolute difference is 
less than 6 percent. 
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These results indicate that, at least for these experiments, the simple scheme de- 
scribed for determining h is adequate for predicting mass ratios. 


SUMMARY OF RESULTS 

A simple model was postulated in NASA Technical Note D-2585 for the problem of 
tank pressurization during outflow. In the present report the equations derived in the 
aforementioned report were further simplified and put into dimensionless form. An 
examination of the dimensionless equations and a large number of numerical solutions 
led to the following results: 

1. The pressurant mass ratio (collapse factor) depends principally on two dimen- 
sionless quantities having the form of modified Stanton numbers, one associated with the 
gas and one with the tank wall. These numbers are functions of tank geometry, outflow 
time, inlet gas properties, tank wall properties, and an average value of heat-transfer 
coefficient. 

2. Other parameters that influence the pressurant mass ratio to a lesser degree are 
the inlet gas temperature, the ratio of liquid surface temperature to inlet gas tempera- 
ture, and the ratio of initial ullage volume to the total volume of liquid outflow. The in- 
fluence of these parameters is negligible in many cases. 

3. A reference Stanton number map is given from which rapid approximate values of 
mass ratio can be found from specified values of reference Stanton numbers. To obtain 
Stanton numbers for a particular problem it is necessary to estimate an average value of 
the heat-transfer coefficient. 

4. For a variety of experimental data the reference Stanton number map was used to 
predict values of mass ratio; average values of heat-transfer coefficient obtained from 
the experiments were used. The predicted values of mass ratio agreed with the experi- 
mental values to within 7 percent, on the average. 

5. Average heat-transfer coefficients were calculated for the experiments by using 
a free convection formula with properties evaluated at the inlet gas temperature. The 
mass ratios predicted using these calculated values of heat-transfer coefficient agreed 
with experimental values to within 6 percent, on the average. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 22, 1964. 
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APPENDIX - SYMBOLS 


c p 

specific heat of gas at constant pres- 

r 

tank radius, ft 


sure, Btu/(lb)(°R) 

r 

dimensionless tank radius, r/(L f - L Q ' 

c w 

specific heat of tank wall, Btu/(lb)(°R) 

st g 

dimensionless gas Stanton number 


integrated average wall specific heat 


(eq. (19)) 


(eq. (23), Btu/(lb)(°R) 

st w 

dimensionless wall Stanton number 

c w 

dimensionless wall specific heat, 
c w /c w 

T 

(eq. (20)) 

temperature of ullage gas, °R 

Gr 

Grashof number 

a 

T 

dimensionless temperature of ullage 

h 

heat -transfer coefficient, 
Btu/(sq ft)(sec)(°R) 

T g 

gas, T/T g 

inlet gas temperature, °R 

L 

height of ullage space, ft 

t l 

liquid surface temperature, °R 

L 

dimensionless height of ullage space, 




L/(L f - L ) 


dimensionless liquid surface tempera- 

L f 

' v f o' 

height of ullage space at end of outflow, 


ture, T L /T g 


ft 

T 

w 

tank wall temperature, °R 

L f 

dimensionless final ullage height, 

T 

dimensionless tank wall temperature, 

L o 

L f /(L f - L q ) 

height of ullage space at start of out- 
flow, ft 

w 

t 

T /T 

w 7 g 

time from start of outflow, sec 

A 

L 

dimensionless initial ullage height, 

t 

dimensionless time from start of out- 

o 

L 0 /(Lf - L 0 ) 


flow, t/tj: 


tank wall thickness, ft 

h 

time at end of outflow, sec 

K 

dimensionless tank wall thickness, 

u 

velocity of ullage gas, ft/sec 


V< L f - L o> 

A 

u 

dimensionless velocity of ullage gas, 

M 

molecular weight 


u/u L 

m 

mass of pres sur ant added during out- 

U L 

velocity of liquid surface, ft/sec 

m i 

P 

flow, lb 

ideal mass of pressurant added 
(eq. (25)), lb 

pressure in tank, lb/sq ft 

X 

X 

space coordinate measured vertically 
downward from top of tank, ft 

dimensionless space coordinate, 
x/(L f - L 0 ) 

Pr 

R 

Prandtl number 
universal gas constant, 

P g 

ullage gas density, lb/cu ft 


ft-lb/(°R)(lb-mole) 

p w 

tank wall density, lb/cu ft 
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